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Abstract: It is shown that perturbations around backgrounds with one non- homogeneous 
dimension, namely of co-homogeneity 1, can be canonically simplified, a property that is 
shown to hold to any order in perturbation theory. Recalling that the problem naturally 
reduces to Id, a procedure is described whereby for each gauge function in Id two Id fields 
are eliminated from the action - one is gauge and can be eliminated without a constraint 
and the other is auxiliary. These results generalize the results of hep-th/0609001 from 
linear to non-linear perturbations and they unify two cases of physical interest: cosmo- 
logical perturbations and perturbations to static spherically symmetric backgrounds. An 
application to black strings is discussed in some detail. 
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1. Introduction 

Analysis of perturbations (gravitational waves) around a given background in General 
Relativity, and in particular making a judicious choice of gauge, is a problem whose com- 
plexity grows with the degree of co-homogeneity of the background, namely the number 
of non-homogeneous coordinates. For co-homogeneity zero, we are dealing with the well- 
understood case of a maximally symmetric space time (flat space, de-Sitter or Anti de- 
Sitter) and once harmonic analysis is employed the problem reduces altogether to algebra. 
The case of co-homogeneity one was solved at the linear level in Q following a concrete 
analysis of the black hole negative mode ^ thereby unifying known results in two cases of 
physical interest: static spherically symmetric space-times (mostly black holes) [||, ||] and 
spatially homogeneous cosmologies || P • The solution for co- homogeneity two and higher 
is unknown and one may expect it to be challenging given the complexity of a specific 
case, that of perturbations around the Kerr spacetime (where Newman-Penrose variables 
are employed). 

Let us briefly review the results of Q. It was shown that once the natural reduction 
to a Id theory is carried out, for each Id gauge function two Id fields can be eliminated: 
one through the gauge and the other being auxiliary. The transformation to the new set 
of fields is local in the non-homogeneous coordinate and invertible. This procedure proves 
that in Id the gauge can be completely eliminated and an invertible transformation exists 
into gauge- invariant fields, which is presumably the reason that the cosmological version of 
this procedure is known as "gauge invariant perturbation theory" . One should stress that 
while it is trivially true that all physical observables including perturbations must be gauge 
invariant it is certainly not true in general that one can transform (one to one and locally) 
into gauge-invariant fields - this property holds only for co-homogeneity one (or zero). 
Actually, one can work in the action formulation and perform a constraint free gauge-fix of 
certain fields (the constraints are contained already in the equations of motion). Moreover, 
not only can the gauge be eliminated thereby eliminating a similar number of fields, but 
also each Id gauge function is responsible for a Id auxiliary field which helps to further 
decouple the action. Accordingly it is said that "the gauge shoots twice". This auxiliary 
field is sometimes considered to be a constraint, though from the current perspective it is 
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certainly not one: a constraint is a residual equation of motion left after fixing the gauge 
while here the gauge gets eliminated altogether. 

The results of H were applied to obtain new insight into gravitational waves in the 
Schwarzschild background (Regge- Wheeler and Zerilli equations) J?]] and to a study of the 
dependence of the black hole negative mode on the space-time dimension (8|. 

One wonders whether the features of linear perturbations carry over to the non-linear 
case and in particular whether each Id gauge function can be continued to be used to 
completely eliminate two Id fields. This question on non-linear perturbations is further 
motivated by applications including cosmology as well as for the computation of the order of 
phase transition associated with the Gregory-Laflamme black string instability || |ll| . 

Experience suggests that the linear part of the equations sets much of their quali- 
tative properties. Indeed, in section ^ we provide and justify a procedure that answers 
affirmatively and to all orders the questions from the preceding paragraph. This work 
can be thought to unify the existing literature on non-linear cosmological perturbations 
|j~2| , [i~3| , |l4|, |l5| with that of non-linear perturbations in spherical symmetry |l6|]. In section 
|3| we spell out in greater detail the application of the general procedure to the non-uniform 
string, but we leave its implementation to future work. 

Comment: This paper is not intended to be submitted to a refereed journal. Papers such as 
this, which suggest a new computational procedure or theory, are easier to referee once they 
are applied and demonstrated in full rather than as an abstract argument, since the former 
confirm the paper's validity and usefulness beyond doubt. Indeed, I plan to apply the 
procedure presented here to the case of the non-uniform black string. On the other hand, 
it is plain that a general procedure together with its rational are self-contained and have 
a value on their own that merits their publication. Thus in order to avoid unnecessary 
difficulties in the process of refereeing I choose to post this paper on the archives, but 
postpone submitting it to a refereed journal until it can be strengthened by an application. 

2. The procedure 

We start this section by stating a canonical procedure to handle non-linear Id actions with 
gauge, and then we continue to explain and justify it. The procedure is 

• Identify the fields and their gauging and thereby the derivatively-gauged (DG) fields, 
namely Id tensors or vectors. 

• It is more economical to fix a gauge before obtaining the action. For each gauge 
function we may eliminate any single field (in whose variation the gauge function 
appears), but not a DG field. This is a constraint-free gauge-fixing Q. 

• Obtain the action (up to the required order in the fields). 

• Write the equations of motion for the DG fields. 

• Use the equations of motion to solve for the DG fields. The solution is obtained 
through series inversion and should be expanded up to the required order. 
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• Substitute the expression for the DG fields back into the action. 

Alternatively, one may skip the gauge-fixing step, and upon substituting the DG 
fields the action will be found to depend only on certain gauge-invariant fields. 

Let us now discuss each item in detail. 
Fields and general set-up. After separation of variables, namely expanding all fields in har- 
monic functions of the homogeneous coordinates, we may reduce the perturbation problem 
to Id by performing the integration over the homogeneous coordinates. More generally we 
consider any Id action with a gauge symmetry. We denote the Id variable 1 by x and the 
fields and gauge fields respectively by the vector notation 

(j) = (j> l (x) i = 1, . . . , uf 

C = C a (x) a = l,...,n G (2.1) 
where np, n G are the number of (real) fields and gauge functions, respectively. 

Identifying the derivatively- gauged (DG) fields. Let us recall the definition of the DG fields 
Q. We assume the gauging to be linear in the gauge functions and at most linear in 
derivatives, which is indeed the case when analyzing perturbations in GR or gauge theory. 
The gauge variation is given by 

6<t> = G 1 t' + G ti (2.2) 

where (Gi(x), Go(x, </>, (j)')) is a pair of np x n G matrices which depend on the coordinate x. 
We allow Go to depend also on the fields <p and their derivatives (non-linear contributions 
to the gauging), but we assume G\ to have no such dependence, which is indeed the case 
in GR perturbations. 

Definition: The image of G\ in field space is called the derivatively- gauged (DG) fields. 

We recognize this algebraic definition to correspond to the standard gauge variation 
of fields which are vectors or tensors from the Id point of view. 

Following the identification of the DG fields, the field space can be split into the DG- 
subspace with coordinates 4>%q and an arbitrarily chosen complement with coordinates 
<j? x , namely 

ft = (0DG> 4>x) a = l,...,n G , r = 1, . . . ,n F - n G . (2.3) 
Obtaining the action. The action 

S = S[4>] (2.4) 

is considered as the input to the procedure. In practice one needs to reduce the perturbation 
problem to Id as already mentioned. To simplify the computation a constraint-free gauge- 
fix should be performed by fixing for each gauge function any field (in whose variation the 
gauge function appears), as long as it is not a DG field. 

In the cosmological set-up x would normally be denoted by t, while in spherically symmetric back- 
grounds it would normally be denoted by r. 
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Writing the equations of motion for the DG fields. This is a straightforward variation of 
the action. 



Solving for the DG-fields order by order and substituting back into the action. The equa- 
tions of motion for the DG fields read 

= -j— = L DG <p DG + C m (c(>x) + O (<f> x , <Pdg) 2 • (2.5) 

OpDG 

This expression highlights the linear part of the equations as obtained in |Q] . We suppressed 
the indices, Log is just a x-dependent no x n G matrix which we assume to be invertible, 
and C m is a linear operator with at most one derivative which mixes (j>x into the equations 
of motion of 4>dg- 

We proceed to perform a change of variables from <Pdg to <Pdg an d substitute it back 
into the action. However, in practice we shall find that it suffices to substitute zero for 
<t>DG thereby somewhat simplifying the procedure. Accordingly we start by explaining the 
change of variables, but for practical purposes one may wish to skip to the paragraph of 

dH). 

We define 

<p DG := LpQ = 4>DG + ^Bg £"m (<t>x) + O ((f>x, <t>DG) 2 (2-6) 



where we used ( |2.5| ) in the second equality to recognize that at the linear level 4>dg ls a 
shifted version of 4>dg- 

In order to substitute the change of variables into the action we need to invert this 
relation and obtain §dg = 4 > DG{ ( t ) DG, ( t ) x)- Since we are dealing with a perturbation 
theory it suffices to perform the inversion perturbatively, and since Ld G is invertible that 
can always be done through a standard procedure as we shall detail shortly. 

Next we substitute the inverted relation into the action. Note that the determinant of 
the transformation is no longer unity, as was the case for linear perturbations, but there is 
no need to account for that as long as classical physics is involved. Actually by construction 
Q2.6|) the equations of motion state that 



WG 



(2.7) 



and hence this sector is guaranteed to decouple, and it suffices to substitute 4>bg 0. 
Note that equations ( |2.7D can be interpreted to define a natural complement in field space 
to the DG fields, rather than the arbitrarily complement X chosen before. Moreover, since 
zero is gauge-invariant so are the 4>dg fields. 

Having realized that eft can be put to zero, let us describe how this could be implemented 
in practice to avoid defining in the first place. We start from the equations of motion 
Q2.5| ), which we now view as an implicit definition for 

4>dg = 4>dg(4>x) ■ (2.8) 

In the context of perturbation theory it suffices to solve for (pne perturbatively. This can 
always be done given that Ld G is invertible (and local in Id) as we proceed to detail. 
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One expands <Pdg as 

4>dg = 4> { S G {4>x) + 4> { dUh) + ■■■ (2.9) 

where the superscript denotes the order with respect to the fields <px- Next one substitutes 
this expansion into the equations ( |2.5[) and solves them order by order. At each order one 
obtains an equation of the form = Ldg4>dg Src^ where the fc'th order source term 

(k) 

depends only upon (fix as well as lower order of 4>dg- Now 4> DG can be solved for by the 
assumed invertibility of Ldg- 

The last step is to substitute ( |2.8| ) into the action. 

Altogether we have arrived at the canonical form of the action and completed our 
description of the process. 



3. Application to the non-uniform string 

In this section we spell out in greater detail how the general procedure applies to the case 
of the non-uniform string. 

The goals for such a computation include 

• The order of the GL phase transition requires expanding the action to fourth order, 
and was preformed in several works (9|, 1C, [llj the last one employing the Landau- 



Ginzburg approach. The current method can be applied to further improve on the 
method by allowing to use a single "master" field, rather than a collection of three 
fields. 

• Performing computations along the non-uniform string branch to a higher order than 
previously achieved. Such computations could be compared to the numerical results 
of JTT| , |l8| . If one is interested only in thermodynamics then it is necessary to compute 
the action up to sixth order to go beyond the existing literature. On the other hand, 
if one is interested in the solutions themselves, only part of the third order was 
computed so far. 

• Improving and extending the results on charged non-uniform strings IS, 2C]. How- 
ever, in this section we shall restrict ourselves to the neutral case for concreteness. 



Order of GL phase transition 

Following is the improved procedure to compute the order of the GL phase transition at 
various space-time dimensions d. 



Set-up and identification of the DG fields. The space time coordinates are the standard 
t, r, z and t runs along the Euclidean time dimension which has length (3 asymptoti- 

cally, r is the radial coordinate, z runs along the compact dimension whose asymptotic size 
is L, and Q-d-2 represents the angular coordinates on the d — 2 sphere. The total space-time 
dimension is d-\- 1. The gauge functions obeying the isometries are £ r , £ z (where £u is the 
parameter of an infinitesimal diffeomorphism) while the DG fields are those which are a 
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vector or tensor from the r coordinate point of view, namely h rr and h rz (where hu U is the 
perturbation to the metric). 

Ansatz. We need an ansatz which keeps the DG fields and it is convenient to take the 
following constraint-free gauge choice 

ds 2 = f{r) dt 2 + fir)- 1 e 2b dr 2 + e 2/3(r) {dz - a dr) 2 + r 2 e 2c dfl 2 d _ 2 , (3.1) 

where 6, a, c are functions of (r, z), (3 is a function of r alone and f(r) = 1 — (ro/r) d ~ 3 . The 
uniform black string is obtained upon setting all the fields to zero. In fixing the gauge we 
note that for z-dependent modes we have two gauge functions, namely £ r , £ 2 which we can 
use to eliminate two fields chosen here to be hu and h zz . However, for the zero mode fields 
the situation is somewhat different Q. As £ 2 and a do not have a zero mode consistent 
with the symmetries only £ r can be used to eliminate a single field, chosen here to be htt 
and hence we retain the zero mode of f3 (or equivalently h zz ). 

Compute the d + 1 dimensional action. Compute the Gibbons Hawking action. This 
could be achieved by computing the Einstein-Hilbert action and integrating by parts to 
avoid terms which include two r-derivatives acting on a single field. For the purpose of 
computing the order of the transition it is enough to retain terms up to fourth order in the 
fields according to the Landau-Ginzburg theory. 

Reduction of the action to Id. Only the z integration is non-trivial and we need to expand 
the fields in harmonic functions. The action has two symmetries which constrain the 
expansion. The first is z — > —z, and the second is z — ► z + L/2 together with e — > —e, 
where e is the small expansion parameter along the non-uniform branch. Both symmetries 
are remnants of the ^-translation symmetry of the uniform string. Expand 

b(r, z) = e b\ cos(kz) + e 2 bo + e 2 62 cos(2kz) 

a(r,z) = ea.\ sin(kz) + e 2 a<i sin(2fc2:) (3-2) 

and c is expanded just like b. All fields are now functions of r alone. Expand the action 
up to the fourth order in e and perform all z integrations. 

Integrate out the DG fields. The DG fields are the harmonics of b, a namely the five Id 
fields 61, 60, 62, «i> «2- For the purpose of this application it is not required to perturbatively 
solve an implicit equation. Next substitute the solutions for the DG fields back into the 
action. Finally, terms with a single derivative can be eliminated in some circumstances 
through integration by parts. 

A Test. The expressions for the quadratic action should agree with those derived in pL 
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(a) (b) 



Figure 1: The expansion of the action up to fourth order in terms of Feynman diagrams. The 
solid four external legs represent the c\ field and the dashed internal propagator in (b) represents 
all possible fields, namely cq, C2, f3. (a) is the quartic vertex of the GL mode while (b) describes the 
contribution from the back-reaction fields. 



namely 

S 2 = S 2 [cx] + S 2 [c 2 ] + S 2 [f3] + S 2 [c ] 

S 2 [ Cl ] = (rf ~ 1} 2 ((i ~ 2) L jdt r d - 2 dr [f c[ 2 - k\\ + V(r) c\] 

= 2(d-l)(d-3)3 

r 2 (2(d - 2)r d - 3 - (d - l)) 2 

S2[P\ = ^L Jdtfr d - 2 dr(3' 2 . (3.3) 

The expression for S 2 [02] is gotten from that of S^fci] by replacing k — > 2k while the 
expression for S 2 [cq] can be deduced from eq. (4.13) in ||. 

First order. Write the equation of motion for c\ and solve the corresponding eigenvalue 
equation to obtain c\ = <t>GL{ r ) and k 2 = \qpy- 

Second order. Solve the back-reaction equations of motion for the fields c 2 , cq, h. 



Substitution in the action. Compute S&(GL) and S 2 (BR) - see [11] for the definitions of 



these notations. Continue to compute C = S<±{GL) — S 2 {BR) which determines the order of 



the phase transition in the canonical ensemble (see also [19|). This computation is encoded 
by the Feynman diagrams of figure |. The results can be compared against table 2 of |lT| . 
It is possible to compute also the coefficient A, allowing full thermodynamic information 
including to pass to the micro-canonical ensemble as in [O]. 
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Higher orders 

The more efficient expressions gotten here should facilitate the computation of higher 
orders. The perturbation theory can be conveniently organized by means of Feynman 
diagrams. For instance, in order to expand the action (or equivalently the free energy and 
obtain the thermodynamics) up to order j we need to compute all tree diagrams (since 
it is a classical field theory) with j external c\ legs, and substitute c\ = (j>GL(r)- In this 
way we compute the free energy F = F{e). Canonically the free energy is a function of 
temperature, or equivalently of k in our case. To obtain F(k) we should compute also 
k = k(e) = fc(°) + e 2 k^ + . . . to the required order. 
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